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We present the first evidence of the decay B^ —>■ t^Vt, using 414 fb^^ of data collected at the 
T(45') resonance with the Belle detector at the KEKB asymmetric-energy e'^e~ collider. Events 
are tagged by fully reconstructing one of the B mesons in hadronic modes. We detect the signal 
with a significance of 3.5 standard deviations including systematics, and measure the branching 



fraction to be B{B 



0.229 



+0.036 



(stat)] 



-^ T'Vr) = (1.79l[5;|^(stat)l!5:5^(syst)) x lO""^. This implies that /s = 
'(syst) GeV and is the first direct measurement of this quantity. 



PACS numbers: 13.20.-v, 13.25.Hv 



In the Standard Model (SM), the purely leptonic de- 
cay B^ -^ T~v^ [l| proceeds via annihilation of h and 
u quarks to a W~ boson (Fig. [T|). It provides a direct 
determination of the product of the B meson decay con- 
stant Jb and the magnitude of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element |Kjb|- The branching 
fraction is given by 



B{B' 



T Vr] = 



"8^ 



1-^ fsWu 



ub\ Tb, 



(1) 

where Gp is the Fermi coupling constant, ttib and rrir 
arc the B and r masses, respectively, and tb is the 
B" lifetime Q]. The expected branching fraction is 
(1.59 ± 0.40) X 10-4 using jKbl = (4.39 ± 0.33) x lO'^, 
determined by inclusive charmless semileptonic B decay 
data 0, TB = 1.643±0.010ps H, and Jb = 0.216±0.022 
GeV obtained from lattice QCD calculations Q . Physics 
beyond the SM, such as supersymmetry or two-Higgs 
doublet models, could modify B{B^ -^ '''"Ri) through 
the introduction of a charged Higgs boson Q. Purely 
leptonic B decays have not been observed before. The 
most stringent upper limit on B~ — > T~i>T comes from 
the BaBar experiment: B{B^ -^ t^Dt) < 2.6 x 10^^ 
(90% C.L.) Q. In this paper, we present the first evi- 
dence for B~ — > T~Dr from the Belle experiment. 

We use a 414 fb~^ data sample containing 449 x 10^ 
B meson pairs collected with the Belle detector at the 
KEKB asymmetric-energy e+e~ (3.5 on 8 GeV) col- 
lider [3 operating at the TUS) resonance {^/s — 10.58 
GeV). The Belle detector Q is a large-solid- angle mag- 
netic spectrometer consisting of a silicon vertex detec- 
tor, a 50-layer central drift chamber (CDC), a system 
of aerogel threshold Cherenkov counters (ACC), time-of- 
flight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented to identify K^ and 



B 



u 




FIG. 1: Purely leptonic B decay proceeds via quark annihi- 
lation into a W boson. 



We use a detailed Monte Carlo (MC) simulation based 
on GEANT Q to determine the signal selection efh- 
ciency and study the background. In order to repro- 
duce effects of beam background, data taken with ran- 
dom triggers for each run period are overlaid on simu- 
lated events. The B~ -^ t~Vt signal decay is generated 
by the EvtGen package [l0|. To model the background 
from e^e" -^ BB and continuum qq {q = u, d, s, c) pro- 
duction processes, large BB and qq MC samples corre- 
sponding to about twice the data sample are used. We 
also use MC samples for rare B decay processes, such 
as charmless hadronic, radiative, electroweak decays and 
b ^ u semileptonic decays. 

We fully reconstruct one of the B mesons in the event, 
referred to hereafter as the tag side (i?tag), and compare 
properties of the remaining particlc(s), referred to as the 
signal side (Bsig), to those expected for signal and back- 
ground. The method allows us to suppress strongly the 
combinatorial background from both BB and continuum 
events. In order to avoid experimental bias, the signal 
region in data is not examined until the event selection 
criteria are finalized. 



The i?tag candidates are reconstructed in the follow- 
ing decay modes: B+ ^ D^*^°n+, W*^^p+, W*^^a+ 
and D'-*^^Ds* ■ The D^ mesons are reconstructed as 



D 



■0 



K 







K^n- 



and K-K+ 



X^ 



XO.o, 



'0, 



r + 



V^boam 



and the D^ mesons are recon- 
structed as D+ -^ KlK+ and K+K--K+. The 7?*° and 
D*^ mesons are reconstructed in D*^ — > £'"7r°,D°7, 
and £)*+ -^ Dfj modes. The selection of Stag can- 
didates is based on the beam-constrained mass Afbc = 
Pg and the energy difference AE = Eb — 
Ehcam- Here, Eb and pB arc the reconstructed energy 
and momentum of the Stag candidate in the e^e~ center- 
of-mass (CM) system, and i?bcam is the beam energy in 
the CM frame. The selection criteria for Stag are defined 
as Afbc > 5.27 GeV/c^ and -80 MeV < Ai; < 60 MeV. 
If an event has multiple i^tag candidates, we choose the 
one having the smallest x^ based on deviations from the 
nominal values of Ai?. the D candidate mass, and the 
D* — D mass difference if applicable. By fitting the Afbc 
distribution to the sum of an empirical parameterization 
of the background shape [11| plus a signal shape p^ . 
we estimate the number of Stag's and their purity in the 
selected region to be 6.80 x 10^ and 0.55, respectively. 

In the events where a -Btag is reconstructed, we search 
for decays of Bgig into a r and a neutrino. Candidate 
events are required to have one or three charged track(s) 
on the signal side with the total charge being oppo- 
site to that of i?tag- The r lepton is identified in the 

e~DpVT. t^~Vt, '!t~ tt'^ Vt and 



five decay modes, /i 



■k~'k^tt~Vt, which taken together correspond to 81% of 
all r decays Q- The muon, electron and charged pion 
candidates are selected based on information from par- 
ticle identification subsystems. The leptons are selected 
with efficiency greater than 90% for both muons and elec- 
trons in the CM momentum region above 1.2 GeV/c, 
and misidcntification rates of less than 0.2%(1.5%) for 
electrons (muons). Kaon candidates are rejected for 
all charged tracks on the signal side. The 7r° candi- 
dates are reconstructed by requiring the invariant mass 



of two 7's to satisfy 
all modes except r" 



M^^ - m^o\ < 20 MeV/c^. For 
— > 7r~7r°i/,-) we reject events with 



it" mesons on the signal side. We place the following 
requirements on the track momentum in the CM frame, 
Pi > 0.3 GeV/cfor/i~z/^^'r ande~z>e^'T,P7r- > 1-OGeV/c 
for 7r~z^T-, Ptt-tt" > 1.2 GeV/c for TT'TT^i/,- andp^-7r+7r- > 
1.8 GcV/c for -k^tt^-k^Vt- We calculate the missing mo- 
mentum of the event in the CM frame (pmiss) from ps 
and the momenta of charged tracks and tt*^ 's on the sig- 
nal side. We require Pmiss > 0.2 GeV/c for ^i^v^Vr and 
e-v^Vr; Pmiss > 1-0 GeV/c for ir^Vr, Pmiss > 1-2 GeV/c 



for n 



and Pmiss > 1-8 GeV/c for tt tt+tt 



In 



order to suppress background where particles produced 
along the beam pipe escape detection, the cosine of the 
angle of the missing momentum (cosSJ^jgg) is required 
to satisfy -0.86 < cos6';;i,3 < 0.95 in the CM frame. 
We further require the invariant mass of the visible de- 
cay products to satisfy |A/^^ — mp| < 0.15 GeV/c^ and 
lAfTTTTTr — TTT-a-l < ^-^ GcV/c^. All thc sclcction criteria 
have been optimized to achieve the highest sensitivity in 
MC. 

The most powerful variable for separating signal and 
background is the remaining energy in the ECL, denoted 
as E'ecl, which is sum of the energies of neutral clus- 
ters that are not associated with either the Stag or the 
7r° candidate from the t~ — *■ 7r~7r''i/,- decay. For neutral 
clusters contributing to i^ECL, we require a minimum en- 
ergy threshold of 50 MeV for thc barrel and 100 (150) 
MeV for thc forward (backward) cndcap ECL. A higher 
threshold is used for the endcap ECL because the effect 
of beam background is more severe. For signal events, 
^ECL must be either zero or a small value arising from 
beam background hits, therefore, signal events peak at 
low E'ecl- On the other hand, background events are 
distributed toward higher i^ECL due to the contribution 
from additional neutral clusters. 

The -Eecl signal region is optimized for each t de- 
cay mode based on the MC simulation, and is defined 
by i^ECL < 0.2 GcV for thc ^~D^i'r, £~i^ei^T and tt";/,- 
modes, and Eecl < 0.3 GeV for thc ■jT~Tr^i'r and 
Tr~ TT~^ Tr~ i^r modes. The -Eecl sideband region is defined 
by 0.4 GeV < Eegl < 1-2 GeV for the /.i^i^^i/,-, e^DgJ^r 
and n~Vr modes, and by 0.45 GeV < E'ecl < 1-2 GeV 
for the TT~Tr^i'T and 7r~7r+7r^z^T- modes. TablcUshows the 
number of events found in the sideband region for data 
(^s°idc) and for the background MC simulation {N^£) 
scaled to the equivalent integrated luminosity in data. 




(GeV) 



FIG. 2: -Eecl distribution for double tagged events, where one 
B is fully reconstructed in the hadronic mode and the other 
B is reconstructed as B~ -^ D*^l~u. The dots indicate the 
data. The solid histogram is the background from BB MC 
{B^ B~ +B'^B'^), while the dashed one shows the contribution 
from B'^B'^ events. 



Their good agreement for each r decay mode indicates 
the validity of the background MC simulation. According 
to the MC simulation, about 95% (5%) of thc background 
events come from BB{qq) processes. Table U also shows 
the number of the background MC events in the signal re- 
gion {N^^). The MC simulation predicts that the back- 
ground in the signal region comes from B~ -^ D^*^^£~D 
semileptonic decays (90%) and rare B decay processes 
(10%). About 30% of the background has K^ candidates 
in the KLM. 

In order to validate the -Eecl simulation, we use a 
control sample of double tagged events, where the -Stag 
is fully reconstructed as described above and Bgig is re- 
constructed in the decay chain, B~ — > D*^l~i) {D*° -^ 
D'^tt'^), followed by 1?° — > K~Tr~^ or K~Tr^TT^TT^ where 
£ is a muon or electron. The sources affecting the E'ecl 
distribution in the control sample are similar to those in 
the signal MC simulation. Fig. [5] shows the -Eecl distri- 
bution in the control sample for data and the scaled MC 
simulation. Their agreement demonstrates the validity 
of the Eecl simulation in the signal MC. 

After finalizing thc signal sclcction criteria, the signal 
region is examined. Fig. [3] shows the i^ECL distribution 
obtained when all r decay modes are combined. One can 
see a significant excess of events in the Eecl signal re- 
gion below Eecl < 0.25 GeV. Table |T] shows the number 
of events observed in the signal region (A^obs) for each t 
decay mode. For the events in the signal region, we ver- 
ify that thc distributions of the event selection variables 
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FIG. 3; -Eecl distributions in the data after all selection cri- 
teria except the one on -Becl. The data and background MC 
samples are represented by the points and the solid histogram, 
respectively. The solid curve shows the result of the fit with 
the sum of the signal (dashed) and background (dotted) con- 
tributions. 



other than -Eecl, such as Mbc and Pmiss, are consistent 
with the sum of the signal and background distributions 
expected from MC. The excess remains after applying a 
K1 veto requirement. 

We obtain the final results by fitting the obtained 
Eecl distributions to the sum of the expected signal 
and background shapes. Probability density functions 
(PDFs) for the signal /s(-Eecl) and for the background 
/bC-E-ECL) are constructed for each r decay mode from 
the MC simulation. The signal PDF is modeled as the 
sum of a Gaussian function, centered at E'ecl = 0, and 
an exponential function. The background PDF, as deter- 
mined from the MC simulation, is parameterized by the 
sum of a Gaussian function and a second-order polyno- 
mial function. The Gaussian function in the background 
PDF addresses deviations from the second-order param- 
eterization, which may arise from a peaking component 
in the lower i^ECL- The PDFs are combined into an ex- 
tended likelihood function, 



£ = 



:,-(nB+nb) 

TV! 



N 



Y[insfsiE,)+nk^iE^)), 



(2) 



where Ei is the -Becl in the ith event, N is the total 
number of events in the data, and Us and n^ are the signal 
yield and background yield to be determined by the fit 
to the whole -Eecl region (0 < E^cl < 1-2)- The results 
are listed in Table HI Table [l] also gives the number of 



background events in the signal region deduced from the 
fit (-/Vb), which is consistent with the expectation from 
the background MC simulation (N^^). 

The branching fractions are calculated as B = Ns/{2 ■ 
£ ■ Nq+b-) where Nq+q- is the number of T{4S) — > 
B'^B~ events, assuming N^+g- — Ng„-ga. The effi- 
ciency is defined as e = e'^'s x e^°^, where e^'^^ is the 
tag reconstruction efficiency for events with B~ — > T~Dr 
decays on the signal side, determined by MC to be 
0.136 ± 0.001(stat)%, and £^°' is the event selection effi- 
ciency listed in Table HI as determined by the ratio of the 
number of events surviving all the selection criteria in- 
cluding the r decay branching fractions to the number of 
fully reconstructed B^ . The branching fraction for each 
T decay mode is consistent within errors. To obtain the 
combined result for all t decay modes, we multiply the 
likelihood functions to produce the combined likelihood 



Ylj^i Cj), and constrain the five signal compo- 



te -^coii 

nents by a single branching fraction. The combined fit 
gives 17.2^47 signal events in the signal region (Ns) and 
24.1^gg in the entire region (tIs). The branching fraction 
is found to be (1.79J:^j^) x 10"^. 

Systematic errors for the measured branching frac- 
tion are associated with the uncertainties in the number 
of B^B^ , signal yields and efficiencies. The systematic 
error due to the uncertainty in N^+g- is 1%. The un- 
certainty in the signal yields arises from uncertainties in 
the signal and background shape, and is determined to 
be ^26 '^^ Here the uncertainty due to the signal shape 
uncertainty is determined by varying the signal PDF pa- 
rameters by the amount of difference of each parameter 
between data and MC for the control sample of double 
tagged events. To determine the background shape un- 
certainty, we vary the Gaussian constant of the back- 
ground PDF by the branching fraction errors from PDG 
for the dominant peaking background sources (such as 
B -^ D^*^°ii^, D° -^ 'k{KYv, etc.). We then add in 
quadrature the variations for the signal and background 
shapes. We take a 10.5% error as the systematic error as- 
sociated with the tag reconstruction efficiency from the 
difference of yields between data and MC for the con- 
trol sample. This value includes the error in the branch- 
ing fraction B{B~ -^ D*'^l~v), which we estimate from 
B{B^ -^ D*^ l^v) in [2] and isospin symmetry. The sys- 
tematic error in the signal efficiencies depends on the t 
decay mode, and arises from the uncertainty in tracking 
efficiency (1 — 3%), 7r° reconstruction efficiency (3%), par- 
ticle identification efficiency (2 — 6%), branching fractions 
of r decays (0.3 - 1.1%), and MC statistics (0.6 - 2%). 
These efficiency errors sum up to 5.6% for the combined 
result after taking into account the correlations between 
the five r decay modes [Tjl • The total fractional system- 
atic uncertainty of the combined measurement is _28%, 
and the branching fraction is 



B{B- 



(1.79+°i^(stat)+°:^6(syst)) X 10- 



The significance is 3.5cr when all r decay modes are 
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TABLE I: The number of observed events in data in the sideband region (A'^s°dc)i number of background MC events in the 
sideband region (A'^'^*^) and the signal region (A'^^'^g ), number of observed events in data in the signal region (A'^obs), number 
of signal (A'^s) and background (A^b) in the signal region determined by the fit, signal selection efficiencies (e°°'), extracted 
branching fraction (JS) for B^ —> t^Vt- The listed errors are statistical only. The last column gives the significance of the 
signal including the systematic uncertainty in the signal yield (S). 



combined, 



where the significance is defined 
where £niax 



as S ~ 
and Lq denote the max- 



-^-21n(£o/'Cmax), 

imum Hkefihood value and likelihood value obtained as- 
suming zero signal events, respectively. Here the likeli- 
hood function from the fit is convolved with a Gaussian 
systematic error function in order to include the system- 
atic uncertainty in the signal yield. 

In conclusion, we have found the first evidence of the 
purely leptonic decay B~ — > t~Dt from a data sam- 
ple of 449 X 10^ BB pairs collected at the T(45) res- 
onance with the Belle experiment. The signal has a 
significance of 3.5 standard deviations. The measured 
branching fraction is (1.79lJ?;^^(stat)to;5i(syst)) x 10"". 
The result is consistent with the SM prediction within 
errors. Using the measured branching fraction and 
known values of Gf, t^b, rrir [3| and tb [3|, the prod- 
uct of the B meson decay constant /s and the magni- 
tude of the CKM matrix element I Vub \ is determined to 



be fB ■ \Vub\ - (10.ll}j(stat)l}j_(syst)) x lO'^ GcV. 



Using the value of lT4b| from [3|, we obtain Jb = 
0.229t^:^i^^(stat)lJ?!|!i^7(syst) GeV, the first direct deter- 
mination of the B meson decay constant. 
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